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1. Introduction 



A systematic way to study bottom baryons from first principles is lattice QCD. Since ami, > 1 
for currently available lattice spacings for large volume simulations, one needs to use for the b 
quark a formalism such as Heavy Quark Effective Theory (HQET). In this paper we consider 
the leading order of HQET, which is the static limit. The spectrum of static-light baryons has 
been studied by lattice methods in the quenched approximation [[lj, ^J, ^] and recently also with 
dynamical sea quarks [Q, H, Bp. Here we use Nf = 2 flavors of dynamical Wilson twisted mass 
quarks. Our intention is to explore the static-light baryon spectrum as fully as possible. Besides the 
experimentally known states A^J, X^/X^, and Q.^ we also predict a number of not yet measured 
states, mainly of negative parity. 



2. Simulation setup 



We use 24 3 x 48 gauge field configurations generated by the European Twisted Mass Collab- 
oration (ETMC). The fermion action is Nf = 2 Wilson twisted mass 1^], ||] tuned to maximal twist, 
where static-light mass differences are automatically &{a) improved JJ). The gauge action is tree- 
level Symanzik improved [10] with j8 = 3.9 corresponding to a lattice spacing a = 0.079(3) fm 



[pp. At the moment we have considered three different values of the twisted mass n q with corre- 
sponding pion masses listed in Table [I]. For details regarding these gauge field configurations we 



refer to [jl2|, |13p. 

In the valence sector we use light quarks, which have the same mass as the sea quarks (corre- 
sponding to u/d quarks), as well as partially quenched quarks with a mass /i q , V ai = 0.0220, which 



is around the mass of the physical s quark Q14[ , |15Q . 





raps in MeV 


number of gauge field configurations 


0.0040 


340(13) 


200 


0.0064 


423(16) 


40 


0.0100 


525(20) 


30 



Table 1: twisted masses /i q , corresponding pion masses mps and number of gauge field configurations 
considered. 



3. Static-light baryon creation operators 

To create static-light baryons, we use operators 

^twisted _ £ afocga^£,(l)^r<^ptwisted^c,(2)^ 

where Q is a static quark operator and %( n > are light quark operators in the so-called twisted basis. 
The upper indices a, b and c are color indices, ^ = 7072 is the charge conjugation matrix and r 1 * 1816 ' 1 
is an appropriately chosen combination of 7 matrices. 
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When discussing quantum numbers of static-light baryons, it is more convenient to first trans- 



form the operators (|3JJ) into the physical basis, in which they have the same structure, i.e. 

£ abcQ U { /^,(l)\r<^pphysical^c,(2) N 



^physical 



(3-2) 



In the continuum the relation between the physical and the twisted basis is given by the twist 
rotation y = exp(iys , v^(o/2)x, where the twist angle ft) = n/2 at maximal twist. At finite lattice 
spacing, however, issues are more complicated: the twist rotation only holds for renormalized 
operators and the QCD symmetries isospin and parity are explicitely broken by 0{a). Nevertheless, 
it is possible to unambiguously interpret states obtained from correlation functions of twisted basis 
operators in terms of QCD quantum numbers. The method has successfully been applied in the 
context of static-light mesons [[if]] and is explained in detail for kaons and D mesons in [17]. For 



details regarding its application to static-light baryons we refer to an upcoming publication [ 18||. 

Since there are no interactions involving the static quark spin, it is appropriate to label static- 
light baryons by angular momentum j and parity of the light degrees of freedom, which are 
determined by rP hysical . Consequently, j = 1 states correspond to degenerate pairs of states with 
total angular momentum 7=1/2 and J = 3/2 respectively. 

To access all possible isospin and strangeness combinations, we consider light quark flavors 

y( 2 ) = ud — du (corresponding to / = 0, S = 0), y/ 1 ) i// 2 ) € {uu ,dd,ud + du) (correspond- 
ing to / = 1, S = 0), y/Wy/ 2 ) e {us, ds} (corresponding to I = 1/2, S = — 1) and y/Wy/ 2 ) = ss 
(corresponding to / = 0, S = —2). 

The operators we have considered and their associated quantum numbers are collected in Ta- 
ble n 



■pphysical 
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/ 
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name 
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name 
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1/2 








K 


1/2 
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YoYs 
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1/2 
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1/2 
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1 





1/2 










1/2 






X 


X 


X 


7o 





1/2 


1 







1/2 









-2 




Yj 


1+ 


1/2, 3/2 


1 







1/2 









-2 


n fe - 


YoYj 


1+ 


1/2, 3/2 


1 





^■b, ^-l 


1/2 









-2 




YjYs 


1 


1/2, 3/2 










1/2 






X 


X 


X 


YoYjYs 


1 


1/2, 3/2 


1 







1/2 









-2 





Table 2: static-light baryon creation operators and their quantum numbers (j'^: angular momentum and 
parity of the light degrees of freedom; J: total angular momentum; /: isospin; S: strangeness); operators 
marked with "X" are identically zero, i.e. do not exist. 



4. Numerical results 

We compute correlation matrices 

Cr, A (0 = (^|<7 ted (0(< isted (0))V)- (4-1) 
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We use several techniques to improve the signal quality, including operator optimization by means 
of APE and Gaussian smearing, the HYP2 static action and stochastic propagators combined with 
timeslice dilution. These techniques are very similar to those used in a recent study of the static- 
light meson spectrum [[| and will be explained in detail in [18]. 

From the correlation matrices (fO|) we extract effective masses by solving a generalized eigen- 
value problem (cf. e.g. j2Q| ] and references therein). Mass values are obtained by fitting constants 
to effective mass plateaus at sufficiently large temporal separations. To exemplify the quality of 
our results, we show effective masses and corresponding mass fits for the A° and the Q.^ baryon in 
Figure [I]. 




2 4 6 8 10 12 
temporal separation 



a 

a 

> 




2 4 6 8 10 12 
temporal separation 



Figure 1: effective masses and corresponding mass fits; left: A° b at fi q = 0.0040 from a 3 x 3 correlation 
matrix; right: Q.^ at /X q = 0.0040 from a 3 x 3 correlation matrix. 



Static-light baryon masses diverge in the continuum limit, because of the infinite self energy of 
the static quark. Therefore, we consider mass differences to another static-light system, the lightest 
static-light meson (quantum numbers = (l/2)~, corresponding to the B and the B* meson, 
which are degenerate in the static limit). In such mass differences the self energies of the static 
quarks exactly cancel. 

We use our results at three different light quark masses (cf. Table [T]) by extrapolating mass 
differences linearly in (raps) 2 to the physical u/d quark mass (raps = 135 MeV). Figure || (left) 
shows this extrapolation for the A°, the X^/X^ and the Q.^ . 

Our results for the static-light baryon spectrum are collected in Table |3} In Wilson twisted 
mass lattice QCD isospin is explicitely broken by 0(a). Therefore, states with 1^0 and different 
I z will have slightly different masses differing by iff {a 2 ). Because of this, in Table ^| two values are 
listed for such states. Note, however, that within statistical errors there is no difference between 
any two such values, which is in agreement with the expectation that isospin breaking effects are 
0(a 2 ) and, hence, should be rather small. 

5. Interpolation in the heavy quark mass and comparison to experimental data 

To make predictions regarding the spectrum of b baryons and to perform a meaningful compar- 
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Figure 2: left: linear extrapolation of static-light mass differences in (mps) 2 to the physical u/d quark mass 
corresponding to raps = 135 MeV; right: linear interpolation between lattice static-light and corresponding 
experimental charm-light mass differences in m c j rag to the physical b quark mass corresponding to m c /mb = 
0.35. 
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name 


m(baryon) — m(B) in MeV 


0+ 










434(46) 
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1 
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1/2 


-1 
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1/2 


-1 




1200(90)/ 1262(77) 
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1/2 


-1 




1233(58)/1285(69) 


1+ 





-2 


a- 


903(42) 


1 
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1315(79) 



Table 3: mass differences of static-light baryons and the lightest static-light meson extrapolated to the 
physical u/d quark mass. 



ison with experimental data, we interpolate linearly m c /m.Q, where mg is the heavy quark mass, to 
the physical point m c /m.f, m(D)/m(B) = 0.35. We do this by using our static-light lattice results 
(cf. Table |3|) corresponding to m c /m.Q = and rather precise experimental data for corresponding 
charm-light systems corresponding to m c /niQ = 1. In Figure || (right) we show the interpolation 
forA^Efr/E^andO^. 

In Table |] we compare results of these interpolations with experimentally available data for b 
baryons. Our lattice results are around 15% larger than the corresponding experimental results, a 
tendency, which has already been observed in our recent study of the static-light meson spectrum 
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[Q, |TJ]. In view of this discrepancy it is interesting to compare with results obtained by other lattice 
groups. When comparing mass differences in units of tq, i.e. dimensionless quantities, e.g. with [Q], 
we find agreement within statistical errors. Note, however, that the scale setting is rather different: 
ro = 0.49 fm in [|J], while the corresponding ETMC value is ro = 0.42 fm. 



name 


m(baryon) — m(B) in MeV (lattice) 


m(baryon) —m(B] 


in MeV (experiment) 




429(30) 
629(28) 
651(28) 
635(25) 
877(27) 


341(2) 
528(3) 
550(3) 
513(3) 
775(7) 


(from [ 
(from [ 
(from [ 
(from [ 
(from [ 


21 
21 
21 
21 

22 


]) 
]) 
]) 
]) 
]) 



Table 4: lattice predictions versus experimental results for mass differences between various b baryons and 
the B meson. 



For the negative parity states listed in Table |3] experimental results for b baryons do not exist. 
There are also no lattice predictions for these states we are aware of. However, one can compare 
with phenomenological models, e.g. with [23], for which one finds qualitative agreement. 



6. Conclusions 



We have studied the static-light baryon spectrum by means of Nf = 2 Wilson twisted mass 
lattice QCD. We have considered three different values of the dynamical quark mass corresponding 
to 340MeV<mps < 525 MeV. We have performed an extrapolation in the light quark mass to the 
physical u/d mass as well as an interpolation in the heavy quark mass to the physical b mass. Our 
results agree within around 15% with currently available experimental results for b baryons. 

Future plans regarding this project include increasing the statistical accuracy of our correlation 
matrices to a level, where also excited states can reliably be extracted. Moreover, we intend to 
investigate the continuum limit, which amounts to considering other finer values of the lattice 
spacing. Finally we plan to perform similar computations on Nf = 2 + 1 + 1 flavor gauge field 
configurations currently produced by ETMC [ 24 ] . 
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